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ABSTRACT 



Aims. In order to search for and study the nature of the low-amplitude and long-periodic radial velocity (RV) variations of massive 
stars, we have been carrying out a precise RV survey for supergiants that lie near or inside the Cepheid instability strip. 
Methods. We have obtained high-resolution spectra of a Per (F5 lb) from November 2005 to September 2011 using the fiber-fed 
Oh . Bohyunsan Observatory Echelle Spectrograph (BOES) at Bohyunsan Optical Astronomy Observatory (BOAO). 

Results. Our measurements reveal that a Per shows a periodic RV variation of 128 days and a semi-amplitude of 70 m s -1 . We find no 
strong correlation between RV variations and bisector velocity span (BVS), but the 128-d peak is indeed present in the BVS variations 
i-C among several other significant peaks in periodogram. 

Conclusions, a Per may have an exoplanet, but the combined data spanning over 20 years seem to suggest that the 128-d RV 
variations have not been stable on long-term scale, which is somewhat difficult to reconcile with the exoplanet explanation. We do 
not exclude the pulsational nature of the 128-d variations in a Per. Although we do not find clear evidence for surface activity or 
rotational modulations by spots, coupled with the fact that the expected rotation period is ~ 130 days the rotational modulation seems 
to be the most likely cause of the RV variations. More observational data and research are needed to clearly determine the origin of 
RV the variations in a Per. 
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^T" \ 1. Introduction investigated. Therefore, finding RV variations in Cepheids-like 

stars near the strip may help to understand these stars. 
^ ! During their evolution, massive stars undergo pulsations by the The bright F supergiant a Per (33 Persei, HD 20902, HR 

*0 . /c-mechanism, which operates in the helium ionization zone. At 1017, SAO 38787, HIP 15863) in the well-studied cluster of the 

O ■ this stage, their expanding outer layers become unstable and they same name has been observed for photometric and spectroscopic 

pulsate regularly within a certain region in the H-R diagram, variations since the end of the 1890s. A history of RV determi- 

. . ■ known as the Cepheid instability strip. Depending on the mass, nations of a Per recorded before the development of precise RV 

stars can move across the instability strip and back again several technique is listed in TableQ] There are few observed RV results 

times as they continue to evolve. Classical Cepheids (Cepheids of a Per despite the long interval. These show that the RV values 

I) and W Virginis stars (Cepheids II) lie at the intersection of vary with a range of ~ 2 km s \ which indicates the possibility 

the strip with the supergiant branch in metal-rich and metal -poor of periodic RV variations. 

stars. In this paper, we present new precise RV measurements of a 

Per acquired in 2005 - 201 1 . In Section 2, we describe the prop- 
Several studies have shown that at least one-half of the stars erties of a Per We ana i yze the RVs and search for variability in 
that reside near the Cepheid instability strip are photometri- Section 3. The origin of the RV variations is discussed in Section 
cally stable at the level of 0.01 ~ 0.03 magnitude (Fernie & 4> and we discuss our results in Section 5. 
Hube 1971; Percy 1975; Fernie 1976; Percy et al. 1979; Percy 
& Welch 1981). However, among the stable non-variable stars, 

low-amplitude variations in the radial velocity (RV) have been 2 . The properties Of the F supergiant a Per 

discovered thanks to high-resolution spectroscopy (Butler 1992; 

Hatzes & Cochran 1995). Despite these efforts, the nature of the The improved investigation of the physical parameters of a Per 

low-amplitude variations and the reason for the existence of non- was reported by Evans et al. (1996). They measured a Per by 

variable stars in the Cepheid instability strip are not yet known, comparing its energy distribution derived from IUE spectra and 

Low-amplitude variations outside the blue edge of the instabil- B, V, R, I, J, and K with model atmospheres convolved with in- 

ity strip may be common, but until now this has been scarcely strumental sensitivity functions. The main parameters were re- 
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Table 1. The RV measurements of a Per recorded before the 
development of high-precision RV measurements. 



RV 


References 


(km s" 1 ) 




-2.40 ± 0.53 


Campbell (1898) 


-2.12 ± 0.4 


Campbell (1901) 


-3.22 ± 0.69 


Vogel (1901) 


-2.69 ± 0.79 


Belopolsky (1904) 


-2.24 


Kustner(1908) 


-1.66 ± 0.61 


Goos (1909) 


-1.57 + 1.43 


Belopolsky (1911) 


-3.42 ± 0.87 


Hnatek (1912) 


-3.0 ± 1.0 


Wilson & Joy (1952) 


-2.4 ± 0.9 


Wilson (1953) 


-2.8 ± 0.5 


Parsons (1983) 


-2.9 ± 0.4 


Beavers & Eitter(1986) 


-2.08 ± 1.62 


Mkrtichian (1990) 



ported by Lyubimkov et al. (2010), who used spectroscopic and 
the photometric methods and improved Hipparcos parallaxes by 
van Leeuwen (2007) to determine the fundamental parameters. 
The compiled basic stellar parameters of a Per are summarized 
in Table [2] 

Stellar radius and diameter are taken from the result of 
Pasinetti-Fracassini et al. (2001). They calculated them with sev- 
eral methods using an intensity interferometer (Hanbury Brown 
et al. 1967) and intrinsic brightness and color (Wesselink 1969; 
Blackwell & Shallis 1977). Because the stellar mass has been 
determined somewhat differently by individual authors, the fi- 
nal mass of a companion (for exoplanet explanation) may vary. 
Table [2] shows the similar determinations of stellar mass even 
though Spaan et al. (1987) estimated a somewhat larger mass of 
8.4 M Q . 

It is difficult to determine a true rotational period because of 
an inaccurate stellar radius and ignorance of the stellar inclina- 
tion. Nonetheless, the values of the projected rotational velocity, 
v rot sin i, in the literature are consistent. Based on the recently 
published value of 20 km s _1 (Lee et al. 2006) and a stellar ra- 
dius of 51 - 55 ^ (Pasinetti-Fracassini et al. 2001), we derived 
a range of the upper limit of the rotational period of 

P rot = 27!7?*/(v rot sin /) = 129- 139 days . (1) 

3. Observations and analysis 

We have acquired 442 spectra of a Per from November 2005 to 
September 201 1 (spanning 133 nights of observations) using the 
fiber-fed high-resolution (R = 90 000) Bohyunsan Observatory 
Echelle Spectrograph (BOES; Kim et al. 2007) attached to the 
1.8-m telescope at Bohyunsan Optical Astronomy Observatory 
(BOAO) in Korea. An iodine absorption cell (I2) was used to 
provide the precise RV measurements. Each estimated signal-to- 
noise ratio (S/N) at the I2 wavelength region is about 250 with 
typical exposure times ranging between 60 and 180 seconds. The 
RV measurements of a Per are listed as online data (Table 4). 

The extraction of normalized 1-D spectra was carried out 
using the IRAF (Tody 1986) software. The I2 analysis and pre- 
cise RV measurements were undertaken using a code called 
RVI2CELL (Han et al. 2007), which was developed at the Korea 
Astronomy & Space Science Institute (KASI). 

Figure [TJ shows RV measurements of a Per and the stan- 
dard RV star r Ceti to demonstrate the long-term stability of the 



Table 2. Stellar parameters of a Per. 



Parameter 


Value 


Reference 


Spectral type 


F5 lb 


Hipparcos 


m v [mag] 


1.82 


Hipparcos 


M r [mag] 


-4.67 


Hipparcos 


B-V[mag] 


0.49 


Mermilliod (1986) 


age [Myr] 


41 


Lyubimkov et al. (2010) 


Distance [pc] 


156 + 4 


Lyubimkov et al. (2010) 


RV [kms- 1 ] 


-2.0 ±0.1 


Gontcharov (2006) 




-3.0± 1.5 


Kudryavtsev et al. (2007) 


Parallax [mas] 


6.44 ±0.17 


van Leeuwen (2007) 




6.43 ±0.17 


Lyubimkov et al. (2010) 


Diameter [mas] 


3.20 ± 0.25 


Koechlin & Rabbia (1985) 




2.8" , 2.9* 


Pasinetti-Fracassini et al. (2001) 


Teff [K] 


6270 ± 120 


Evans et al. (1996) 




6240 ± 20 


Lee et al. (2006) 




6541 ± 84 


Kovtyukh (2007) 




6350 ± 100 


Lyubimkov et al. (2010) 


[Fe/H] 


- 0.28 ± 0.06 


Lee et al. (2006) 




-0.07 ±0.09 


Lyubimkov et al. (2010) 


l°g£ 


1.5 


Evans et al. (1996) 




0.58 ± 0.04 


Lee et al. (2006) 




2.0 


Kovtyukh et al. (2008) 




1.90 ± 0.04 


Lyubimkov et al. (2010) 


R* [Re] 


60.7 


Spaan et al. (1987) 




55°, 51* 


Pasinetti-Fracassini et al. (2001) 


M* [M ] 


7.2 


Parsons & Bouw (1971) 




7.2 


Becker et al. (1977) 




8.4 


Spaan et al. (1987) 




7.3 ± 0.3 


Lyubimkov et al. (2010) 


U [L Q ] 


5500 


Spaan et al. (1987) 


v rot sin i [km s~'] 


17.9 ± 1.0 


Gray & Toner (1987) 




18 


Hatzes & Cochran (1995) 




17.1 ± 1.0 


de Medeiros et al. (2002) 




20 


Lee et al. (2006) 


P rot / sin ! [days] 


129-139 


Derived 


Vmicro [km s -1 ] 


3.20 ± 0.05 


Lee et al. (2006) 




5.3 ± 0.5 


Lyubimkov et al. (2010) 


Vmacro [km S" 1 ] 


10 


Hatzes & Cochran (1995) 



Intensity interferometer ( *' Intrinsic brightness and color 



BOES. The figure shows that the RV of t Ceti is constant with an 
rms scatter of 6.7 m s _1 over the timespan of our observations. 

We calculated the Lomb-Scargle periodogram of the BOAO 
RV time series of a Per, which is a useful tool to investigate long- 
period variations for unequally spaced data (Lomb 1976; Scargle 
1982). Figure [2] shows a significant power at f\ = 0.0078 cd _1 
(P = 128.2 days). We find a false-alarm probability (FAP) of < 
10~ 6 by a bootstrap randomization process (Kiirster et al. 1999). 
We could find neither periods of 87.7 days (Hatzes & Cochran 
1995) nor 77.7 days (Butler 1998). 

4. The origin of the RV variations 

Giants and supergiants exhibit pulsations as well as surface ac- 
tivities, resulting in low-amplitude RV variability on different 
time scales. Short-term (hours to days) RV variations have been 
known to be the result of stellar pulsations (Hatzes & Cochran 
1998), whereas long-term (hundreds of days) RV variations with 
a low-amplitude may be caused by three kinds of phenomena: 
stellar oscillations, rotational modulations, or planetary com- 
panions. To establish the origin of the low-amplitude and long- 
periodic RV variations of a Per we examined the Ca II H & K 
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Fig. 1. RV measurements of the RV standard star r Ceti (top panel) and a Per from November 2005 to September 201 1 (bottom 
panel). 
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Fig. 2. Lomb-Scargle periodogram of the RV measurements of a 
Per. The periodogram shows a significant power at a frequency 
of 0.0078 c d _I corresponding to a period of 128 days (top panel) 
after subtracting the main frequency variations (bottom panel). 
Arrows indicate previously reported periods of 87.8 and 77.7 
days. 
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Fig. 3. The spectral line of the Ca II H & K region in a Per. 
It displays three spectra exhibiting positive, zero, and negative 
RV deviations, respectively (top panel), and a combined figure 
of the three features (bottom panel). No emission features are 
seen around the Ca II H & K central region, and any systematic 
differences in line profiles are invisible. 



lines, the Hipparcos photometry, spectral line bisectors, and fit- 
ted a Keplerian orbit to RV data. 



4.1. Ca II H& Klines 

We used the averaged Ca II H & K line profiles at the upper, near 
zero, and lower part of the RV curve to check for any systematic 
difference in the Ca II H & K region related to the RV varia- 
tions. We selected spectra spanning a long time interval: the first 
positive deviation on JD-245375 1.125798 (15 January 2006, RV 
= 144.2 m s" 1 ), near zero deviation on JD-2454535.956330 (10 
March 2008, RV = - 7.0 m s _1 ), and the last negative deviation 
on JD-2455495. 383283 (25 October 2010, RV = - 109.3 m s" 1 ) 
of RV curve. Figure |3]shows that the Ca II H & K cores show no 
emissions at the line centers, and there is also no systematic dif- 
ference among the line profiles at the three extreme parts of the 
RV curve. It means that F5 lb a Per has exhibited no detectable 
peculiarities in Ca II H & K lines relevant to RV variations. This 
agrees with Pasquini et al. (2000), who found that supergiants in 



their sample of stars show no signatures of Ca II H & K chromo- 
spheric activity. 

4.2. Hipparcos photometry 

We analyzed the Hipparcos photometry of a Per to search for 
possible brightness variations relevant to the period of 128 days. 
For three years, the Hipparcos satellite obtained 87 photometric 
measurements of a Per. Figure |4] shows the Lomb-Scargle pe- 
riodogram of the measurements. There are no significant peaks 
near the frequency of f\ = 0.0078 c d -1 corresponding to the pe- 
riod of 128 days. The Hipparcos data maintained a photometric 
stability down to rms of 0.0048 magnitude (0.26%) during those 
three years. 

4.3. Line bisector variations 

The analysis of a line shape can help to determine the origin of 
the observed RV variations (Queloz et al. 2001). The RV varia- 
tions by planetary companions should not produce any changes 



3 



Byeong-Cheol Lee et al.: Detection of the 128 day radial velocity variations in the supergiant a Persei 



Q_ 




0.02 0.03 

Frequency (eld) 



Fig. 4. Lomb-Scargle periodogram of the Hipparcos photomet- 
ric measurements of a Per for three years. The arrow marks the 
location of the period of 128 days. 



in the spectral line shape, whereas the surface inhomogeneities 
or pulsations do. Stellar rotational modulations of surface fea- 
tures or non-radial pulsations can create variable asymmetries in 
the spectral line profiles. Thus, the accurately measured varia- 
tions in the shapes of stellar lines may provide valuable infor- 
mation about surface inhomogeneities, pulsations, or planetary 
companion. The bisector velocity span (BVS) is defined as the 
difference between two bisectors at the top and bottom of the 
profile and is used to quantify the changes in the line profile. 

We measured the BVS using the least-squares deconvolu- 
tion (LSD) technique (Donati et al. 1997; Reiners & Royer 
2004; Glazunova et al. 2008). We used the Vienna Atomic Line 
Database (VALD; Piskunov et al. 1995) to prepare the list of 
spectral lines. A total of ~ 3100 lines within the wavelength re- 
gion of 4500 - 4900 and 6450 - 6840 A were used to construct 
the LSD profile, which excluded spectral regions around the I2 
absorption region, hydrogen lines, and regions with strong con- 
tamination by terrestrial atmospheric lines. 

We estimated the BVS of the mean profile between two dif- 
ferent flux levels of central depth levels, 0.8 and 0.25, as the span 
points. The BVS as a function of JD is shown in Figure [5] We 
calculated the Lomb-Scargle periodogram of the BVS shown in 
Figure [6] We see a significant peak frequency 0.0078 c d~' (a 
FAP ~ 0.5%) relevant to the period of 128 days. There are in- 
deed several significant peaks with higher or comparable power 
in the periodogram. Although it is difficult to claim with cer- 
tainty whether the signal near the 128-d period is physically real, 
it is unlikely that noise would cause this peak. The variations in 
the BVS might have a higher magnitude compared with the sig- 
nal at 128 days and hinder it. 



4.4. Keplerian orbital fit 

We analyzed our BOAO data using the Lomb-Scagle peri- 
odogram and found the best-fitting Keplerian orbit with a P — 
128.2 + 0.1 days, a semi-amplitude K - 70.8 + 1.5 m s , and 
an eccentricity e — 0.10 + 0.04. The RV curve phased to the 
period of 128 days is shown in Figure [7] and it exhibits clear 
variations. We derive the minimum mass of a companion m sin i 
= 6.6 + 0.2 Mj up at a distance of a = 0.97 AU from a Per, only 
four times the stellar radius R+ = 51 - 55 R Q (~ 0.25 AU). 

As can be seen in Figure [7] the rms of the RV residuals are 
44.8 m s _1 , which is larger than the typical internal error of in- 
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Fig. 5. BVS variations of 
September 201 1. 
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Fig. 6. Lomb-Scargle periodogram of BVS variations of a Per. 
The arrow marks the location of the period of 128 days 

dividual a Per RV accuracy (~ 20 m s~'). We also noticed some 
systematic pattern in the periodogoram of RV residuals after re- 
moving the best fit. Figure [2] (bottom panel) shows the Lomb- 
Scargle periodogram. We found a FAP of ~ 10~ 3 at the high- 
est peak of the residual RV measurements, but there are several 
peaks at this level, which suggests noise. Irregular periodic RV 
variations in a Per appear with a 2K amplitude of ~ 150 m s _I in 
a short period of time. The scatter seems to be caused by varia- 
tions due to the intrinsic stellar variability or stellar oscillations. 
All the orbital elements from the BOAO data set are listed in 
Table[3] 



5. Discussion 

F-G yellow supergiants are very important for studying the 
Galactic structure owing to their high luminosity and young 
age (Klochkova & Panchuk 1988; Giridhar et al. 1998). 
Furthermore, all types of supergiants are good science cases in 
pulsations regarding the Cepheid instability strip. According to 
the the traditional criteria (Iben & Tuggle 1975; Fernie 1990), 
a Per is located outside the blue edge of the Cepheid instability 
strip (Hatzes & Cochran 1995; Bulter 1998). 

If supergiants with luminosities and temperatures similar to 
the Cepheids do not pulsate in the Cepheid instability strip, 
they may be related to non-variable supergiants (NVSs), which 
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Fig. 7. RV measurements of a Per phased to the orbital period of 
128 days obtained at BOAO. The solid line is the orbital solution 
that fits the data with an rms of 44.8 m s _1 . 

Table 3. Best-fitting Keplerian orbit of the hypothetical planet a 
Perb. 



Parameter 


Value 


Period [days] 


128.2 ±0.1 


K[m s _1 ] 


70.8 ± 1.5 


^periastron [JD] 


2451513.00 ±4.91 (1999.912) 


e 


0.10 ±0.04 


(D [deg] 


208.37 ± 12.36 


/(«) [M e ] 


(4.6473) x 10-" 


tr(O-C) [m s- 1 ] 


44.8 


with M = 7.3 [M Q ] 




m sin i [M Sup ] 


6.6 ± 0.2 


a [AU] 


0.97 



were investigated by Fernie & Hube (1971). The presence of 
the NVSs in the instability strip is quite enigmatic. s-Cepheids 
are pulsating stars with low amplitude variations of the light, 
color, and RV. They can be considered as objects in the stage 
between normal Cepheids and the NVSs in the instability strip. 
Andrievsky & Kovtyukh (1996) found no significant differences 
in the C and O abundances between them and concluded that 
these two groups of stars have undergone the first stage of core 
helium burning and are at the same evolutionary stage. From the 
photometrist's point of view, a Per can be classified as an NVS 
near or on the blue border of the Cepheid instability strip. 

In this work, we found compelling evidence for a low- 
amplitude and long-periodic RV variation in the F supergiant 
a Per. A strongly periodic variability in RV measurements was 
detected with a period of 128 days. We will examine possible 
origins of these RV variations in the next subsections. 



5.1. Pulsations 

The estimated period of the fundamental radial mode pulsations 
is several days. Therefore, the 128-d period of periodic varia- 
tions found in a Per cannot be caused by classical /(-mechanism 
pulsations operating in Cepheids. On the other hand, we found 
additional 119-d RV variations on top of the low-amplitude 
Cepheid variations of 3.97-d period in Polaris (Lee et al. 2008), 
a Cepheid that lies inside the instability strip. Lee at al. (2008) 
speculated about the existence of an yet unrecognized mecha- 
nism that drives low-amplitude and long-period oscillatory vari- 
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Fig. 8. Lomb-Scargle periodogram of the data combined with 
the McDonald (Hatzes & Cochran 1995) and Lick observatory 
(Butler 1998) RV measurements of a Per. Arrows indicate ~ 128 
and ~ 77 days. 
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Fig. 9. Lomb-Scargle periodograms of the BOAO, McDonald, 
and Lick observatory RV measurements of a Per. The peri- 
odograms show a frequency of 0.0078 c d~' at the BOAO data 
set (solid line) and 0.0077 c d" 1 at the BOAO + McDonald + 
Lick observatory data set (dashed line), respectively. 



ations of the radius in F-supergiants and in all classical Cepheid 
stars. Lee et al. (2008) expressed the need for a detailed com- 
parative study of the precise RV for a sample of F-type super- 
giant stars including Cepheids that can provide additional con- 
straints on the nature of the long-term radius variations detected 
in Polaris, a Per is the second supergiant star in our survey that 
shows nearly identical (in terms of RV period and amplitude) 
long-period quasi-periodic variations and seems to support the 
assumptions made by Lee et al. (2008) about the common na- 
ture of the oscillatory variations of the radius in F-supergiants. 



5.2. Rotational modulations 

Hatzes & Cochran (1995) and Butler (1998) argued that the 
cause of the RV variations of a Per is not likely to be the or- 
bital motion, but rather the intrinsic stellar variations such as ro- 
tational modulation or non-radial pulsations. The 119-d period 
with a 2K amplitude of ~ 280 m s _1 variations in the residual 
RV measurements in the F supergiant Polaris could also have 
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Fig. 10. RV measurements of a Per phased to the orbital period of 129.9 days. Filled circles, triangles, and open circles denote the 
RV measurements of BOAO, McDonald, and Lick observatory, respectively. The solid line is the orbital solution that fits the data 
with an rms of 58.5 m s . 



been caused by rotational modulation of surface features (Lee et 
al. 2008), for example. 

Supergiants typically have rotational periods of several hun- 
dred days (Rao et al. 1993). a Per shows a projected rotational 
velocity of 17-20 km s _1 (see the Table 2) corresponding to 
a period of 129 - 139 days (upper limit), very close to the ob- 
served period of 128 days. It is also close to the mean projected 
rotational velocity of 15 km s _1 for 16 F supergiants (Danzinger 
& Faber 1972). Although the real rotational period is likely to 
be shorter by the factor of sin i than 129 - 139 days, it does not 
eliminate the rotational origin of RV variations. When the incli- 
nation of a Per is very close to 90 degrees there is a chance that 
the observed RV variations are caused by the rotational modula- 
tion. 

a Per is an object with weak longitudinal magnetic fields 
of Bi = 0.82 + 0.37 G (Grunhut et al. 2010) but it exhibits the 
strongest and most complex Stokes V profile of all 30 super- 
giant stars observed. Such a low value of magnetic fields does 
not mean, however, that the local fields are weak as well. The 
local magnetic fields can be stronger and have a complex sur- 
face structure. Even a weak magnetic field in spots is enough 
to suppress the convection motion and make the velocity field 
(macroturbulence) inside the spot lower than the surrounding 
photosphere. This would alter the line shapes and result in an 
RV shift. As shown for Polaris by Hatzes & Cochran (2000) 
such macroturbulent spots can account for RV variations of up 
to ~ 100 m s _1 . Such spots would produce no detectable chro- 
mospheric activity indicator (Figure|3} or photometric variations 
(Figure 0). We need a direct confirmation of the existence of the 
local magnetic fields and the surface inhomogeneities in a Per. 

5.3. Planetary companion 

The important features of the hypothetical planetary RV varia- 
tions are the strong period, shape, and phase stability of varia- 
tions. We clearly found 128-d variations. However, the measure- 
ments in the previous records by Hatzes & Cochran (1995) and 
by Butler (1998) did not show the same variations. Nevertheless, 
the data combined with previous records do not show a signifi- 
cant power, just two somewhat high peaks exhibited near ~ 128 



and ~ 77 days, with a moderate power of ~ 20 (Figure|8]l. If we 
have regular, long-term stable 128-d variations in two data sets, 
the combination of them should increase the power in the result- 
ing periodogram, which is not the case. We also calculated the 
Lomb-Scargle periodogram of the whole RV measurements of 
BOAO, McDonald (Hatzes & Cochran 1995), and Lick observa- 
tory (Butler 1998). Figure |9]shows the periodogram of the com- 
bined data against that of BOAO. We find a significant power at 
at /[ = 0.0077 cd 1 (P = 129.9 days) from the whole data set, 
approximately the same frequency in BOAO data, and an FAP 
of this periodicity is less than 10 . However, the Lomb-Scargle 
power decreases, which means that 128-d variations have not 
been stable on the long-term scale. This is a strong criterion to 
reject the exoplanet hypothesis. The RV measurements of the 
whole data set phased to the period of 129.9 days are shown in 
Figure [10] It does not conform perfectly with whole RV mea- 
surements. 

Assuming an exoplanet (least possible) hypothesis, we can 
estimate the minimum mass for the planetary companion of 6.6 
Mj np with an orbital semi-major axis of 0.97 AU and an eccen- 
tricity of 0. 1 . 
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17 


.7 


4216.945450 


-90.5 


26 


.3 


5580. 


.917689 


-46.7 


21. 


.7 
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Table 4. continued. 



JD-2 450 000 


ARV 




JD-2 450 000 


ARV 


±cr 


JD-2 450 000 


ARV 




day 


m s 


m s ^ 


day 


m s ^ 


m s ^ 


day 


m s 


m s 


-i 


3980. 


.240219 


47.9 


17.1 


4216.948054 


-85.8 


23.8 


5580.918858 


-45.0 


1 o 

iy 


o 

a 


3980. 


.244617 


20.9 


20.6 


4216.950658 


-42.5 


24.9 


5580.920027 


-48.2 


1 "7 


A 

y 


3981. 


.219074 


11.7 


21.3 


4262.310982 


99.4 


18.1 


5580.921196 


-42.3 


O £ 

Z3. 


-7 


3981 


.222825 


31.8 


21.9 


4262.312915 


153.0 


25.3 


5580.922365 


-53.3 


1 


n 

. 1 


3985. 


.309197 


-10.8 


17.3 


4263.309016 


88.1 


25.2 


5581.060109 


-76.7 


1 ft 


.1 


3985. 


.312863 


11.9 


17.0 


4263.311204 


131.3 


29.3 


5581.061730 


-75.7 




o 

.1 


4007. 


.097216 


12.4 


20.2 


4264.305744 


73.0 


40.9 


5581.063246 


-99.8 


Zl. 


.0 


4007. 


.102367 


44.1 


30.5 


4268.307138 


146.1 


37.0 


5581.064762 


-60.7 


10. 


A 

.4 


4008. 


.283631 


32.8 


20.5 


4268.310032 


163.2 


23.5 


5581.066278 


-78.8 


1 "7 


.0 


4008. 


.289400 


40.4 


19.6 


4269.305285 


110.8 


19.2 


5583.991513 


-78.8 


1 n 
19 


o 

.1 


4008. 


.300638 


41.6 


15.8 


4269.306766 


95.2 


26.6 


5583.992681 


-47.7 


1 Q 


-7 


4008. 


.304955 


39.3 


14.1 


4274.301293 


166.1 


24.4 


5583.993850 


-52.1 


1 1 . 


A 

y 


4008. 


.309176 


20.7 


19.8 


4274.303562 


147.3 


20.8 


5583.995019 


-43.9 


la. 


A 

.9 


4011 


.172797 


26.2 


17.5 


4274 305807 


175.5 


41.0 


5583 996200 


-54.1 


OA 

10. 


A 
.0 


4011 


.176356 


8.8 


19.4 


4382.068443 


64.9 


22.2 


5583.997369 


-49.4 


i n 


A 

u 


4011. 


.180045 


48.5 


19.1 


4382.073177 


65.5 


21.8 


5583.998537 


-65.2 


i n 

iy 


A 

.9 


4012. 


.258615 


41.9 


13.9 


4382.077969 


52.3 


22.9 


5583.999706 


-47.6 


i n 

iy 


.4 


4012. 


.262224 


38.9 


20.2 


4382.176494 


50.0 


21.5 


5584.000875 


-59.3 


1 a. 


.4 


4012. 


.265720 


37.0 


15.0 


4382.181112 


45.5 


24.5 


5584.002044 


-43.4 


O 1 

Zl 


o 
.5 


4012. 


.269076 


29.9 


17.0 


4388.088347 


36.0 


21.4 


5585.911355 


-53.3 


OQ 

15. 


A 

.4 


4014. 


.353056 


46.6 


19.4 


4388.166801 


38.0 


21.6 


5585.912524 


-37.8 


oo 

11. 


A 

.9 


4014. 


.358959 


76.1 


23.7 


4388 272849 


44.3 


20.8 


5585 913693 


-53.1 


oo 

11. 


.5 


4014. 


.370534 


54.0 


16.2 


4389.279420 


59.6 


17.3 


5585.914862 


-61.9 


O A 

24, 


t 

,1 


4014. 


.380673 


112.1 


16.9 


4389.356812 


50.5 


17.7 


5585.916031 


-50.6 


O 1 

11. 


/r 
.0 


4015. 


.156985 


64.7 


16.8 


4389.996033 


36.9 


18.7 


5585.917200 


-42.7 


OA 

10. 


.3 


4015. 


.160457 


61.3 


17.0 


4390.047540 


59.6 


20.4 


5585.918368 


-49.5 


OA 

10. 


A 

.4 


4015. 


.163987 


51.6 


20.0 


4390.120449 


39.2 


20.3 


5585.919537 


-69.1 


i n 

iy 


A 

.4 


4015. 


.167333 


62.7 


15.9 


4390 156099 


26.1 


19.4 


5585.920706 


-86.8 


OA 

10 


.3 


4018. 


.042674 


117.5 


22.6 


4390.199516 


53.7 


20.0 


5585.921875 


-56.8 


OA 

10 


.3 


4018. 


.048010 


99.1 


21.8 


4390.225096 


50.8 


19.7 


5588.904606 


-57.6 


1 o 

la 


o 

.1 


4018. 


.053774 


69.9 


20.0 


4390.263304 


28.4 


22.5 


5588.906122 


-57.1 


1 o 

la 


.0 


4018. 


059161 


58.0 


19.3 


4390 285636 


36.6 


18.8 


5588 907638 


-47.0 


1 H 

1 /. 


-7 


4018. 


.337553 


67.6 


17.7 


4390.310292 


26.0 


22.2 


5588.909154 


-55.8 


1 O 

lo. 


o 

.1 


4018. 


.340979 


73.6 


16.6 


4390 333389 


46.7 


19.5 


5588.910671 


-40.2 


1 n 

iy. 


A 

.9 


4018. 


.344405 


54.9 


17.9 


4390.356098 


51.4 


14.9 


5588.912430 


-47.9 


1 "7 


.1 


4023. 


.111979 


80.9 


16.4 


4393.106831 


71.3 


17.5 


5588.913946 


-37.0 


1 a 


o 

.a 


4023. 


115690 


81.1 


16.0 


4396.305248 


76.5 


19.4 


5588.915462 


-58.9 


1 H 


A 
.0 


4023. 


.119370 


67.1 


18.3 


4396.309137 


97.6 


17.9 


5588.916978 


-56.4 


OA 

10. 


A 

.4 


4023. 


.122982 


62.2 


15.6 


4420.134262 


111.9 


27.2 


5588.918494 


-52.2 


1 a 


.0 


4027. 


.084998 


44.2 


25.5 


4420.136947 


79.7 


25.2 


5590.082499 


-62.0 


OA 

10. 


.5 


4027. 


.089269 


63.0 


14.4 


4420.139870 


93.7 


27.4 


5590.084362 


-57.2 


O A 

Z4. 


Q 

.a 


4027. 


.092950 


61.1 


19.8 


4424.088364 


53.6 


26.9 


5590 086226 


-48.3 


oo 

11. 


A 

.9 


4027. 


.096804 


63.1 


20.5 


4424.090181 


35.1 


24.1 


5590.088089 


-83.1 


o c 

ZD. 


A 
.0 


4027. 


.101180 


64.8 


20.9 


4453 078318 


-52.7 


21.4 


5590.089952 


-67.7 


oi 
15 


A 

u 


4036. 


.149659 


6.3 


20.7 


4453.081524 


-59.9 


22.8 


5590.091815 


-66.3 


O 1 

Zl. 


A 

.4 


4036. 


.153158 


20.3 


16.7 


4453.084730 


-71.7 


19.8 


5590.093679 


-36.7 


0^2 

15. 


A 
.0 


4036. 


.156353 


20.4 


17.3 


4454.893158 


-71.3 


22.8 


5590.095542 


-71.7 


OA 

ZD. 


A 

y 


4036. 


.159513 


17.6 


16.0 


4454.895669 


-45.8 


23.4 


5590.097405 


-54.1 


oo 

11. 


Q 

.a 


4038. 


.221880 


26.0 


21.1 


4454.898181 


-67.8 


20.5 


5590.099269 


-53.2 


O^l 

15 


Q 

.a 


4038. 


.225627 


7.8 


21.2 


4467.937579 


-141.5 


19.3 


5614.912587 


-59.5 


in 
5 1 . 


.1 


4038. 


.229181 


4.5 


16.5 


4467.939454 


-130.6 


18.8 


5614.914115 


-84.8 


0^ 

ZD. 


A 

u 


4043. 


.185848 


-8.6 


16.3 


4472.022170 


-93.3 


23.5 


5614.915978 


-69.2 


1 / 


A 

.9 


4043. 


.189115 


-12.5 


14.6 


4472.026533 


-93.1 


20.4 


5614.917841 


-66.0 


OO 

11. 


o 
.1 


4043. 


.192437 


-23.4 


20.9 


4472.030954 


-103.8 


19.4 


5614.919716 


-65.4 


1 A 

iy. 


A 

y 


4047. 


.304018 


-13.0 


23.6 


4472.034981 


-86.4 


22.0 


5614.921626 


-61.4 


1 "7 

1 /. 


.3 


4047. 


.307501 


-15.7 


17.6 


4483.086998 


-35.6 


15.2 


5614.923489 


-68.8 


1 "7 
1 /. 


.3 


At\^ 1 
4U31. 


.UooOlZ 


-30.6 


21.4 


4483.089266 


-20.8 


14.0 


5614.925352 


-71.1 


i n 

iy 


A 

.9 


4051. 


.095221 


-34.5 


24.6 


4483.091153 


-21.0 


18.3 


5614.927227 


-71.4 


27 


.3 


4051. 


.102073 


-77.2 


21.2 


4504.962102 


117.0 


23.9 


5614.929090 


-56.9 


19 


.9 


4052. 


.089730 


-56.8 


19.9 


4505.091755 


89.6 


19.1 


5819.167964 


-4.7 


19. 


8 


4052. 


.093607 


-55.3 


19.5 


4505.094093 


107.9 


18.6 


5819.169168 


2.8 


19. 


6 


4052. 


.097554 


-53.0 


24.2 


4505.952087 


88.9 


19.2 


5819.170360 


6.0 


15. 


.5 


4071 


.069923 


-114.5 


23.5 


4505.954796 


101.2 


20.0 


5819.171564 


13.8 


18. 


.3 


4071. 


.071752 


-116.6 


21.2 


4505.957307 


106.6 


20.7 


5819.172757 


-8.4 


20. 


8 


4071. 


.073940 


-108.4 


18.7 


4520.923815 


6.3 


24.2 


5819.173949 


10.2 


22. 


.8 


4071. 


.076428 


-106.7 


21.9 


4520.926350 


9.9 


22.0 


5819.175164 


-13.5 


19 


.6 
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Table 4. continued. 



JD-2 450 000 


ARV 




JD-2 450 000 


ARV 


±cr 


JD-2 450 000 


ARV 


±<x 


day 


m s 


m s _1 


day 


m s _1 


m s _1 


day 


m s 


m s _1 


4071.078963 


-121.2 


19.1 


4535.956330 


-7.0 


18.2 


5819.176414 


-1.2 


24.2 


4081.276441 


-59.0 


17.1 


4535.958309 


-33.0 


22.8 


5819.177618 


-4.5 


20.9 


4081.280226 


-66.7 


20.4 


4537.991880 


-21.7 


30.5 


5819.178822 


3.4 


17.8 


4081.283525 


-67.0 


21.5 


4537.993894 


-22.7 


29.7 


5820.164120 


15.8 


22.5 


4082.041085 


-54.5 


24.9 


4545.940856 


-14.5 


34.7 


5820.165637 


-3.7 


19.1 


4082.044233 


-51.4 


22.1 


4546.925839 


3.6 


31.3 


5820.166609 


-24.2 


20.6 


4082.047833 


-59.9 


24.0 


4546.928350 


5.0 


30.4 


5820.167570 


-0.5 


18.4 


4090.035309 


-45.6 


26.4 


4552.955671 


24.7 


14.0 


5820.168542 


-10.7 


16.3 


4090 038908 


-12.6 


18.1 


4552.957581 


6.2 


20.0 


5820.169514 


-18.7 


20.4 


4092.996622 


2.8 


18.7 


4552.959502 


13.9 


17.0 


5820.170487 


-15.1 


18.7 


4093.000684 


10.5 


20.4 


4640.287985 


120.3 


26.1 


5820.171418 


4.2 


22.4 


4098.901678 


9.6 


24.1 


4640.290566 


105.3 


20.8 


5820.172333 


-0.1 


19.7 


4098.903757 


3.5 


19.4 


4640.293148 


109.5 


24.1 


5820.173247 


-3.0 


16.3 


4098.905398 


12.5 


20.2 


4640.295763 


115.0 


23.7 


5820.174162 


-14.4 


19.6 


4098.906847 


-3.4 


22.6 


4717.135611 


-103.7 


25.4 








4104.131330 


44.3 


26.5 


4717.138721 


-117.8 


24.9 









